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Abstract 

We demonstrate that the -matrices in the superintegrable chiral Potts model possess the 
Onsager algebra symmetry for their degenerate eigenvalues. The Fabricius-McCoy comparison 
of functional relations of the eight- vertex model for roots of unity and the superintegrable chiral 
Potts model has been carefully analyzed by identifying equivalent terms in the corresponding 
equations, by which we extract the conjectured relation of Q-operators and all fusion matrices 
in the eight-vertex model corresponding to the TT-relation in the chiral Potts model. 
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Introduction 

The aim of this paper is to show the Onsager algebra symmetry of r^-matrices of the chiral 
Potts model (CMP) [14] in the superintegrable case. Inspired by the discussion in [29] on the TQ- 

(2) 

relation of eight-vertex model, we observed in [39] that the degeneracy of Tp -eigenvalues occurs 
at the superintegrable point p. Indeed in [29] , Fabricius and McCoy introduced the fusion matrices 
of eight- vertex model for the "root of unity" 

?j=^, (gcd(m 1; iV) = l), (1) 
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and proposed the conjecture on functional relations related to the Q-matrices, of which the struc- 
tures strikingly resemble the set of functional equations revealed in the study of TV-state CPM in 
[14]. In this paper, we make a detailed analysis on the comparison of these two theories (superin- 
tegrable CPM and eight-vertex model for roots of unity) through the equivalence of operators and 
quantities appeared in functional equations in an "identical" manner (for the correspondence, see 
table (39) of this article). The effort enables us to further explore one more conjecture on the rela- 
tionship between Q-operators and the whole set of fusion matrices, the equation (40) of this paper, 
corresponding to the TT-relation in CPM which has played the fundamental role in derivation of 
functional equations in [14]. In the special case where the elliptic nome vanishes, the eight- vertex 
model has the six- vertex limit for the root of unity fj = in which case there is the s/2-loop 

algebra symmetry [20, 25, 26, 27]. Hence the degeneracy of the six- vertex transfer matrix is mainly 
contributed from finite-dimensional representations of s^-loop algebra, which are determined by 
evaluation parameters as the roots of the Drinfeld polynomial. However the essence of Drinfeld 
polynomial is encoded in the Q-operator, by the analysis made in [28] on the eight-vertex Q-matrix 
in connection with the function in [21] on degeneracy of the transfer matrix. The analogy between 
eight-vertex model and CPM, though only loosely related in the present stage, still enables us to 
suggest that the degeneracy of r^-operator in the superintegrable CPM, previously indicated in 
[39], should further posses a certain symmetry structure similar to the s^-loop algebra symmetry 
in six-vertex model in the root of unity case. In this paper we show that it is indeed the case. 
We identify the full symmetry algebra for t^, hence all matrices to be the Onsager algebra, 
a renowned Lie algebra appeared in the seminal paper of Onsager on the free energy solution of 
the two-dimensional Ising model [36], (for the Onsager algebra, see [18, 19, 23, 38] and references 
therein). In the study of (finite-dimensional) representations of Onsager algebra, a proper and use- 
ful realization of Onsager algebra is to identify it with the Lie-subalgebra of s^-loop algebra fixed 
by a standard involution, by which the s^-loop algebra representations (with certain constraint on 
evaluation parameters) give rise to all the irreducible representations of Onsager algebra [18, 38]. 
From the Onsager algebra structure hidden in the Hamiltonion derived from superintegrable CPM, 
we identify its role in connection to the symmetry of t^ 2 ) -matrices. Furthermore, through the func- 
tional relations and eigenvalue spectrum of the transfer matrix in the superintegrable CPM, Onsager 
algebra indeed describes the full symmetry for the degeneracy of t( j ) -matrices, where only spin-| 
representations occur in the related s^-loop algebra representations. Hence, the Onsager algebra 
symmetry of -matrices observed in this article, in contrast with the s/2-loop algebra symmetry 
of six-vertex model, further enhances the common features shared by both superintegrable CPM 
and eight-vertex model for the roots of unity cases. 

This paper is organized as follows. In Sec. 1, we briefly review the main features of the iV-state 
CPM (e.g., [3, 5] and references therein), summarize the set of functional equations of the model 
in [14], then describe the specific form of r^ 2 ) -matrix in [39] and formulae in the superintegrable 
case considered in this paper. In Sec. 2, we begin with the formalism and conjectured functional 
relations in eight- vertex model of roots of unity in [29] , then make comparison of all entities appeared 
in functional relations of the superintegrable CPM and eight-vertex model in a precise and one- 
to-one corresponding manner, by which a conjectured functional equation in eight-vertex model 
corresponding to TT-relation in CPM naturally ascends. The efforts encourage us to propose 
a speculation about functional relations in a hypothetical solvable lattice model with the similar 
structure as in the superintegrable CPM. Based on the common features appeared in superintegrable 
CPM and eight-vertex model for roots of unity, we outline the procedure and "ingredients" for a 
such speculated general scheme. The next two sections consist of main results of this paper on the 
symmetry algebra of superintegrable CPM. In Sec. 3 we start with some basic facts in Onsager 
algebra and its relation with superintegrable CPM and the quantum Z^-spin chain Hamiltonian 
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in [32, 34]. Using the explicit forms of monodromy matrix in the definition of and the Onsager 
algebra generators in the Z^-spin chain Hamiltonian, we show the Onsager algebra symmetry 
of r^-matrices in sup erintegr able CPM. While in Sec. 4 we discuss the Bethe equation of t^- 
matrix and the explicit results on r^) -eigenvalues to which the Onsager algebra symmetry of 
t( j ) -matrices can be understood from the Bethe Ansatz approach based on eigenvalue spectrum of 
the superintegrable CPM transfer matrix. Using the explicit form of eigenvalues of CPM transfer 
matrix, hence the energy expression of the Z^r-spin chain Hamiltonian, one can determine that only 
spin-| representations occur in the Onsager algebra symmetry of t( j ) -matrices with the explicit 
T ti) -eigenvalue expression. The mathematical nature of P-polynomial associated to the Bethe 
equation of the theory has been rigorously investigated so that the full symmetry structure for the 
degeneracy of t^) -matrices is revealed in Onsager algebra representations. We close in Sec. 5 with 
some concluding remarks. 

Notations. To present our work, we prepare some notations. In this paper, Z,R, C will 
denote the ring of integers, real, complex numbers respectively, Z^r = Z/iVZ, and i = For 
N > 2, we fix the N th root of unity, 

2iri 

uj = e N , 

and C N is the vector space consisting of all TV-cyclic vectors with the basis {| n )} ng Zjv - Define 
X, Z the operators of by the relations, X\n) = \n + 1), Z\n) = uj n \n) for n 6 Z^r, then the 
Weyl relation and iV th -power identity property hold: XZ = lj~ 1 ZX, X n = Z N = 1. For a positive 

integers n, we denote by <§> the tensor product of n-copies of the vector space C N . 



1 The iV-state Chiral Potts Model 

We start with the summary of functional equations of chiral Potts model in [14], (for a short 
version of explanations, see [39]), then conduct the discussion in the superintegrable case for later 
consideration. This will also serve to establish the notation. 

In the study of two-dimensional iV-state CMP, the "rapidity" of the statistical model is described 
by a four-vector ratios [a, b, c, d] in the projective 3-space P 3 satisfying the following equivalent sets 
of equations, 

/ ka N + k'c N = d N , f a N + k'b N = kd N , 

W: \ kb N + k'd N = c N . ^ \ k'a N + b N = kc N , (2) 

where k,k' are parameters with k 2 + k' 2 = 1, and k' ^ ±1,0. The above relations define 2U as 
an algebraic curve of genus N 3 - 2N 2 + 1. We shall confine our discussion of CPM only on the 
full homogeneous lattice by taking p = p' in [14]. Hereafter, elements in P 3 will be denoted by 
p,q,r, - ■ ■ etc., and we shall use the variables x, y, fj,, t, A to denote the following component-ratios 
for an element [a, b, c, d] G P 3 , 

a b d ab d N N 

x:=-, y:=-, fi := -, t := - (= xy) , A := ^ (= M ). 

Associated to an element p in P 3 , the coordinates and variables will be written by a p , b p , x p , t p ■ ■ ■ 
so on whenever if it will be necessary to specify the element p. The variables, x, y, //, form a system 
of affine coordinates of P 3 , by which an equivalent form of the defining equation for the curve 2H 
in (2) is 

kx N = 1 - k'n~ N , ky N = 1 - k'n N , (x, y, fi) G C 3 . (3) 
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The rapidity curve 2U has a large symmetry group (with the order 4N 3 for N > 3 [4, 39]), in which 
the following automorphisms will be needed in later discussions, 



U : [a, b, c, d] i-> [ua, b, c,d], ((x, y, p) i-> (lux, y, fi)j , 
C : [a, b, c, d] i-> [b, a, d, c], y(x, y, fi) h-> (y, x, fi' 1 )^ . 



(4) 



Note that U N = C 2 = 1. By eliminating the variable fj, N in (3) , 2U becomes a iV-fold unramified 
cover of the genus (N — l) 2 curve, 

x N + y N = k(l + x N y N ) . (5) 

By (3), the variables, (t, A) = (t p ,X p ) for p € 2U, define the following hyperelliptic curve of genus 
iV-1, 

WV : t = p , (6) 

which is a iV-ramified quotient of the curve (5). The variables t, A in the above curve have played a 
major role of rapidities for solutions of various problems in CPM, see e.g. [9, 13, 35] and references 
therein. 

The Boltzmann weights W Piq , W Ptq of the TV-state CPM are defined by coordinates of p, q 6 2U 
with the expressions: 

W Pl q(n) _ i-jn d v b q -a v c q u)i I _ , fj, p -, n yr n y q -u)ix p 
W p , q (0) llj'=l b p d q -c p a q u>3 I V?' Vp-w'x q 

W p , q {n) _ ytn ua v d q -d v a q ui I _ , \ n T\ n ux v —u)ix q 

W p , q (0) ~ ll 3 =1 c p b q -bpc q ^ y ~ WpVq) llj=l ^-o^j/p 

By the rapidity constraint (2), the above Boltzmann weights have the iV-periodicity property for n. 
Without loss of generality, we may assume W Ptq (0) = W Ptq (0) = 1. On a lattice of the horizontal size 
L with periodic boundary condition, the combined weights of intersection between two consecutive 

rows give rise to an operator of <£> C N , which defines the transfer matrix of the TV-state CPM: 

L 

T p (q) a , a > = HWp^cn - a\)W p ^{ai - a' l+1 ) , p,q£W, (7) 
1=1 

where a = (a±, . . . , &£) > 0"' = ■ ■ ■ , &l) with ai,a[ G Z^r. By the star-triangle relation of 
Boltzmann weights: 

N-l 

J2 W v{o' - n)W pr (a - n)W pq (n - a") = R pqr W pq (a - a')W pr (a' - a")W qr (a - a") 

n=0 



with R pqr = ^jrf- and f pq = y ^^-i^ ^ J > one has the commutativity of transfer matrices 
for a fixed p G 22J: 

[T p (q), T p (q')}=0 , q, q' G 2XJ . 

By (3), T p (q) for a fixed p depends only on the values of (x q ,y q ), parametrized by the curve (5). 
Hence we shall also write T p (q) by T p (x q ,y q ) whenever this is convenient. The Z^r-spin operators 

X, Z of give rise to the system of Weyl operators of ® C^: 

X j = l®---®X jth ®l®---®l, Zj = 1® ■■■®Z jth ® 1® (l<j<L) 
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with ZiXj = uj^XjZi, [Zi,Zj] = [Xi,Xj] = and Zf = Xf = 1. The spin-shift operator of 

<8> C N will be again denoted by X = Ylf =1 Xj if no confusion could arise. Then T p (q) commutes 
with X and the spatial translation operator Sr, which takes the j'th column to (j + l)th one for 
1 < j < L with the identification L + 1 = 1. The eigenvalues of X and Sr, denoted by e 1( ^, e lP for 
Q £ Zjy, P £ Zl, are quantum numbers, called the Z^-charge and total momentum respectively. 
We denote 

f p (q) = T p (q)S R . 

In the study of CPM as a descendent of the six-vertex model, Bazhanov and Stroganov dis- 
covered a five-parameter family of solutions of Yang-Baxter equation associated to the six-vertex 
R-matrix with entries in terms of operators X, Z acting on "quantum space" [15, 31]. In 

particular, there is a Yang-Baxter solution G(t) associated to an element p = [a, b, c, d] £ P 3 with 
the form [39]: 

2 _ / 6 2 - td 2 X (be - ujadX)Z \ 

b G p (t) - ^ _ t{bc _ adx)z -! _ tc 2 + J , t€C, (8) 

which satisfy the Yang-Baxter relation: 

R(t/t')(G p (t)(g)l)(l<g)G p (t')) = (l®G p (0)(Gp(t) <g>l)R(t/t') , (9) 

aux aux aux aux 

where R(t) is the following matrix of 2-tensor of the "auxiliary space" C 2 , 

( tuj - 1 \ 

t- 1 u-1 

t(u-\) w(t-l) 

\ toj - 1 J 

By the auxiliary-space matrix-product and quantum-space tensor-product for a finite size L, the 
operator 

L 

(g) G p/ (t) = G p ,i(t) ® ■ ■ ■ ® G P)L (t) , := G p (t) at site £ , (10) 

again satisfies the Yang-Baxter relation (9), hence the traces, tr aux 0^ =1 G p ^(t) (t £ C), form a 

L pj (2) 

family of commuting operators of <8> C . The t p -operator is defined by 

L 

T ( p 2) (t)=tr aux (®G p/ (ujt)) for tGC. (11) 

(2) 

When p £ 2U, the transfer matrix T p (q) of CPM can be derived from t p (t q ) as the auxiliary 
"Q" -operator in a similar way as in the TQ-relation discussion of eight- vertex model in [7]. One 
has the following r^T-relation ((4.20) in [14], or (14) in [12]) 1 : 

T^{t q )T p {ujx q ,y q ) = <p p (q)T p (x q ,y q )+Tp p (Uq)XT p (uj 2 x a ,y q ) , (12) 

where <p p (q) := ( (yp gg ) iy ) ) £ , W P (<?) := ( ^g^r^ ^' and U is the automOT P hism of W 
in (4). An equivalent form of the r^T-relation (12) is the following expression of T p 2 ^ in terms of 

tfKh) = (v P {Q) T p{xq,yq) +^p(Uq)XT p (uj 2 x q ,y q )^T p (uJx q ,y q y 1 , (13) 



R(t) 



1 Tp 2 \q),T p (q) in this paper are the operators r^ g , T 9 in [14] respectively. 
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by which, the commutativity of T p (g)'s ensures 

[4 2 \t q ),T p (x ql ,y ql )] = , for p,q,q' £ 2U . (14) 

There are operators r (j) (t) for < j < N with r^ 0) (t) := 0, T^\t) := I, which are related to 
transfer matrices T p (q) by the following TT-relations ((3.46) for (l,k) = (j,0) in [14], or (13) in 
[12]): 

T p (x q ,y q )T p (y q ,Jx q ) = r p , q h,, p , q ^\t q ) + *('«)*MJ " ^N-j)^ (15) 

p n ^ • ^ /\r v, - N ( x p~ x <i)(yp~y<i)( t p ^tg ^ L i, _ ( rii -1 Up( x p- U ' mx i) \ L j 

tOTU - J - JV ' Wllere r P>9 - I (^-^)(^-< )(tp _ ig) I , %;p,g - I 1 lm=l (y p -u>™x q ){t p -u™t q ) J and 



*w = ( ^ ( ;-r t)2 ) L (=^(^ P w) 



M\) - (^(i-fc'A^ J (- ^ ^^•..■ ; y-,-; ) 
In particular, the relation (15) for j = N (or 0) becomes 

'Ny$ N -*(y p -y q )(y p -x q )^ 



(16) 



rp i yf, ,_( ^Vp \Vp-Vq)KVp-Xq) Y {N)(+S n7 A 

J-pKXq.yqjlpKyq^Xq) - y ( y N _ yN^yN _ X N^ J T P V r <?J- \ LI ) 



!p Uq J\Up ~q 

By (14) and (15 ), the following "fusion relations" hold for r^'s ( (4.27) of [14] ): 

T$\t)T?\tj-H) = zi^-^Xr^it) + 4 J+1 \t), 1 < j < N , 
T ^ +1 )(t) : = z^Irf" 1 ^^) + u(t)J , 



(18) 



where u(t) = a p (X) + a p (A _1 ), and z(t) in (16). Therefore the fusion operators t p ^ can be con- 
structed recursively from T p 2 ^ and relations in (18) by setting t p °^ = 0, = /, and one can 
express T p \t) as a "polynomial" of of degree (j — 1) with coefficients in powers of X. Indeed 
the explicit form of T p \t) for 2 < j < N + 1 in terms of T p 2 ^ and X is as follows [39]: 

r<%)=nVv*)+x;'(-*) 1 e (n ,S r > ){u " t] )^ ] 

s=0 k=l l<u<'i 2 <'-<'i fc <j-2 V £=l *)7p s=0 7 

where ^ <' means i£ + 1 < v+i- 

Using (13) and (18), one can successively express t p ^ in terms of T p (q), then obtain the following 
r^T-relations for 1 < j < N + 1, ((4.34) in [14]): 

4 j \q) = T p (x q ,y q )T p (^x q ,y q )Y:^lo {^ P {q)^p{Uq) ■ ■ ■ ^ p {U m - 1 q)x 

(20) 

p p (U m+1 q) ■ ■■^ p (W- 1 q)T p (u; m Xq,yqr 1 T p (u; m + 1 Xq, yq )- 1 X3- m - 1 ' 

Indeed, by II^Jo 1 fpi^l) = a p(\) an d Y\f=oWp{U^q) = a p {\ q l ), r^T-relation with the fusion 
relation is equivalent to the r^T-relations : (12) + (18) <=> (20). By (17) and (20) j=N , one 
obtains the functional equation of chiral Potts transfer matrices T p ((4.40) of [14]): 

fp(yq,x q ) = E™=oOn ;P (g)T p (x 9 ,y 9 )T p ^ (21) 



6 



where C m , p {q) = ^ P {q)^ P {Uq) ■ ■ ■ M^' 1 q)%(U m+1 q) ■ ■ ■ y P (^ JV ~ 1 g)( ^^)(g^f 1 ) L - Note 

that substitutions of (20) in (15) all give the same relation (21). Hence under the scheme of t^T- 
relation (12), the ^-functional equation (21) is equivalent to TT-relations (15). The functional 
equations (18) and (15) have provided an effective method to study the eigenvalue problem of CPM 
by reducing the relations in forms with the variables (t, A) in (6) using the following formula ((2.40) 
in [14]): 

Hereafter, the chiral Potts model considered in this paper will only be the superintegrable case, 
i.e., the vertical rapidity p is given by 

i 1 — k' ± 

P ■ x p = y p = rj2 , n p = 1 , where r\ = fc/ ) n . (23) 

We shall omit the subscript p of the operators t p \ T p (x q ,y q ), T p ( Xq ,q q ) for p in (23), e.g. 
means Tp \ etc. Then monodromy matrix in (8) becomes 

r(f) _ ( G° (t) Gj(t) \ _ / 1- V -HX r,^(l -uX)Z \ , . 

and r^ 2 ) -matrix (11) is expressed by 

r (2) (t) =t w ((g)G,)M) = ( £ (Gi)S(G 2 )^---(G L )2i)(^) (25) 

i=\ ^ ai,a2,...,ai ' 

where = 0, 1 for all The fusion relation (18) becomes 



T^{t)T^\^-H) = (1 - ^uj'J-H) 21 T^~ l \t) uj l X + r^' +1 )(t), 1 < j < TV 

T (N+i) {t) = (1 _ r/ -i^2L T (iV-i)(^) LJ Lx + 2(1 - rT^) L . 



(26) 



By (19), one has 

T », t )=n^(^) + E(-^ e (n j;j;;;;l? ( ll t) n^(^))P7) 

o— n i w»- / 'Ia.sa n ' "V' V"^ ^ / o— n / 



s=0 fc=l l<ii<'i 2 <'-<'ifc<j-2 x £ : 



for 2 < j < TV + 1. 
Define 



T( Xq ,y q )(l-rj^ X ^ T( Xq ,y q )(l- V ^ X ^ 

(4c P {q) - — — — , ^ C pW) - — — — • (28) 

N L (1 - r]— Xq ) L N L (l - r]— Xq ) L 

Since li^flM = T i^y q ) ^L x yi by ^ Q cp ( q )N and Q cp ( g )JV depend on i y on the value 

(1-77 "3" 0JX q ) L (l-V^~Xq) L 

of (t q ,X q ) in the curve (6). Hence, up to a iVth root of unity, we may write Q C p(q), Qcp(q) simply 
by Q C p(t q , X q ), Qcp(tq,X q ) as well. Then r^T-relation (12) becomes 

r (2 \t q )Q cp {LJt q , X q ) = (1 - r]-Hq) L q cp {t q , X q ) + (1 - ^Lot^Q^uhq, X q )uj L X , (29) 
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and r^T-relations (20) have the form 

m=o ^ 'it 



q cp {ujn q ,\ q )q cp (u m+ H q ,\ q )-\u L xy- m - i y (30) 



The TT-relations, (15) and (17), become 

-n_ (H^y / i_ 



Qc P (^,A,)Q cp (a;^ 9 ,A- 1 ) = V" J ( n^i"^)^ ^ 

+ !E^ ^ ) ) ' (31) 

Q cp (t g , A,)Q cp (t ? , A" 1 ) = ^(i 5 ). (32) 

By (21), one has the following functional equation of Q C p{o)' 



(I _ n~ N t N ) L 7V_1 i L n _1_m 

Qc P (Cg) = 1 /V"£ E n _^-U )L Qep( g )Qc P (^ g )- 1 Qep(^ +1 g ) 

iV m=0 ^ W ^ V 



(33) 



where C, U are the transformations in (4) . 

2 Eight-Vertex Model for the Root of Unity Cases 

In [29], Fabricious and McCoy observed the similar structures between the eight- vertex model for 
the "root of unity" cases and the superintegrable TV-state CPM. They derived and proposed a set 
of functional equations for the eight-vertex model analogous to those of CPM in [14] as follows. 

The transfer matrix for the eight- vertex model in the root of unity r/ in (1) with even L columns 
and the periodic boundary condition is 2 

T (2) (v)^ u = TrWsQu^WsijMi, ^) ■ ■ ■ W % (ji L , v l ) 

where the Boltzmann weights are parametrized by elliptic theta functions 3 , 

W 8 (M)|i,i =W 8 (1,1)| M =@(2?j)e(v-rj)H(v + 7j) =a(v), 

W 8 (-W)|i,i =W S (1, l)|_i,_i =Q(2rj)H(v-rfie(v + rfi =b(v), 

W 8 (-l,l)|i,-i =^ 8 (l,-l)|-i,i =H(2rj)@(v-rj)@(v + rj) =c(v), 

W 8 (l,-l)|i,-i =W 8 (-1,1)|_ 1>1 =H(2ri)H(v-rj)H(v + ri) =d(v). 

In [6], Baxter defined the "auxiliary" Q-operators, Qr{v) and its companion Ql{v), to study the 
eigenvalue problem of (v) through the following TQ-relations: 

T^(v)Q R (v) = h L (v - fj)Q R {v + 2rf) + h L (v + ?i)Q R (v - 2rj), 
Ql(v)tW(v) = h L {v- v)Ql(v + 2rj) + h L (v + ^)Ql(v - 2rj), 



2 The L, N in [29] are changed to N, L for the parallel notations of the chiral Potts model used in this paper. 
3 Here we use the standard conventions for Jacobi theta functions H(v), Q(v) with elliptic integrals K, K' of nome 

= e^: H(v) =2Y:» 1 (-irV"- 1/2)2 sin &£&L, Q(v) = 1 + 2 ^~ 1 {^Y^ «* z™. 
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where h(v) = <d(0)@(v)H(v). Furthermore, the following relations hold: 

T(v + Nrj) = T(v), h(v + Nrj) = h(v), 

Ql(v)Qr(v') = Q R (v')Q L (v), Q r , l (v + Nrj) = S m ^Q R , L (v) = Q R , L {v)S^ 

where S = ]Jj =1 aj. Define Qsv( v ) = Qr^Qr^o)^ 1 = Ql(^o) _1( 5l('w) for a base element vo. 
Then {Qsv(v)}v and {T^ 2 \v')} v i are commutative families of operators such that [T^ 2 \v'), Qsv(v)] = 
and the following TQ-relation holds: 

T^^Qsviv) = h L {v - rj)Q 8V (v + 2rj) + h L {v + rj)Q 8V {v - 2rj) , 

which is equivalent to 

T^(v + 2?j)Q R (v + 2rj) = h L (v + Srj)Q R (v) + h L (v + rj)Q R {v + 4fj) . (34) 

In [29], Fabricius and McCoy defined the "fusion matrices" T^\v) recursively by setting T^\v) = 
0, T^(v + 2rj) = h L (v + rj), and the following fusion relations ((3.4) (3.5) (3.6) and (3.15) in [29]) 4 : 

T®(v + 2?j)T^(v + 2jrj) = h L (v + (2j + \)if)T^-^(v + 2rj) + T^ +1 \v + 2rj), l<j<N, 
T( N+1 \v + 2rj) = T( N -^(v + 4rj) + 2h L (v + rj)S mi 

(35) 

which, by (34), are compatible with the T^Q-relations ((3.7) in [29]): 
j-i 

T^\v + 2rj)=J2 hL ( v + ( 2m + ^)Qr{v)Q~ r {v + 2m?))Q R (v + 2jrj)Q R 1 (v + 2(m + l)ff) . (36) 

m=0 

Furthermore, they proposed the conjectured Q-functional equation ((3.2) in [29]): 

JV-l 

Q L (v-iK') = e^rQ L (v )AQ R {v ) ]T h(v + (2m + l)rj) L Q R 1 (v + 2m7j)Q R (v)Q R 1 (v + 2(m + l)?i) , 

m=0 

(37) 

which is equivalent to the following one ((3.10) in [29]): 

Q R {v)[e^Q L {v )AQ R {v )^ 1 Q L (v - \K') = T^ N \v + 2^)S m ^ (38) 

Then the (r 2 T, fusion, r^T, (TT) N , Q cp -functional) relations (29) (26) (30) (32) and (33) in CPM 
can be matched to the corresponding ones in eight- vertex model, (34), (35), (36), (38) and (37), 
respectively in an identical manner by paring the quantities in the following table: 





Chiral Potts model 


Eight vertex model 


u 




v i — * v -\- 2rj 


V 


C 


v ^ v — \K' 


t( 2 )(w) 


r( 2 )(t) ^ 


(v+2rj) 
h L (v+rj)h L (v+3rj) 


Q(w) 


Qc P (q) 


Qr(v) 


ip(w) 


(l-ri-%) 


h(v + T])^ 1 


t^(w) 




T<J)(v+2rj) 

ni-i ^ L (^+(2fc-i)?) 


Q(w) 


Qc P (q) 




A 


cj l X 


J 


B 






Mi 


(i-r,- N t») L 


7TiL( t , + 2jl ? ) , , , 

e 2K Q L (v )AQ R {v ) 





(39) 



4 For the comparison with functional relations in chiral Potts model, the T"'- fusion relations and Q_R-functional 
equation of eight- vertex model in [29] are changed to the equivalent forms here by using the transformation, v i— > v+2rj. 
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Furthermore, by (39) one arrives the following conjectured relation in eight-vertex model corre- 
sponding to the set of TT-relations (31) in CPM: 

Q R {v) (e" iL ™ 2iv) Q L (vo)AQ R (vo?) ^ Q L (v - iK' + 2jrj) = T&(y + 2rj) + T^\v + 2(j + l)?j)S mi 

(40) 

for < j < N. By using (36), each relation in (40) is equivalent to the functional equation (37) as 
in the CPM case. 

By the comparison made above on superintegrable CPM and eight-vertex model for roots of 
unity, the discussion suggests that the set of functional equations might possibly occur as well for 
some general models in the following manner, by which when applying to CPM and eight-vertex 
model, the procedure produces all the functional equations in previous discussions. Consider a 
solvable lattice model with Boltzmann weights in a spectral curve W, which is 2B in (2) (or Wk> 

-ttK' 

in (6)) in CPM , and the elliptic curve with nome q = e K in eight- vertex model. The theory 

is built-up from families of operators of "quantum space" ® C m , {t( 2 )(w), Q(w), Q(w)} we vK> such 
that Q(w) and Q(w) are non-degenerated with Q(w)Q(w') = Q(w')Q(w), and the following tQ- 
relations hold: 

t( 2 )(w)Q(Uw) = ^(w) L Q(w) +^(Uw) L Q(U 2 w),4 , 

Q(Uw)t( 2 )(w) = <p(w) L Q(w) + ^(Uw) L Q(U 2 w)A [ ' 

Here (p is a rational function on W, U is an automorphism of W with the following " iV-periodic" 
property with respect to t^ 2 ) and Q, Q: 

t^fU^w) = t^(w), Q(U Jv w) = Q(w)B, Q(XJ N w) = Q(w)B, 

and A, B are operators of <S> C m , commuting with all Q(w) and Q(w), with [A, B] = 0, A N = 
B 2 = 1. Define the operators t^(w) (2 < j < N + 1) by the t^Q-relation: 

t 0-)( w ) = £ ^ !B^^^ Q( w )Q(um w) -lQ (U j w) Q (ir +l w) -l A j-m-lY (42) 

Then t^'s satisfy the fusion relations by setting = and = 1: 

t^)(w)t( 2 )(U^- 1 w) = ¥ >(U^- 1 w) 2L t^- 1 )^) A + t^' +1 )(w), 1 < j < N , 

t ( JV + 1 )(w) = <p(w) 2L t^-^CUw)^ +2(^nf=o^(U fc w)) B . (43) 

The fusion matrices t^'s are expected to relate to Q-operator by the following QQ-relation, 

Q(w)Mr 1 Q(VU i w) = - i 1 -tk')(w) + 7^-^ -t^-^tPw) , < j < N, (44) 

L 

where V is an order 2 automorphism of W, and Mj are operators of ® C m , such that by using 
(42), all relations in (45) reduce to a single Q-functional equation: 

N-l A-l-m 

Q(Vw) = M ^ (TJmw)L Q(U m w)- 1 Q(w)Q(U m+1 w)- 1 . (45) 

m=0 ^ \ ' 

The equivalent notions with those in CPM and eight vertex model are given in table (39). Then the 
functional relations among Q and Q give rise to the corresponding ones in these two models. 
Note that among all these relations, tQ-relation (41) and QQ-relation (45) are the basic ones for 
the theory as indicated in the study of the CPM, however is still unknown for QQ-relation in 
eight-vertex model case. 
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3 Onsager Algebra Symmetry of r^-matrices 

The Onsager algebra is the infinite-dimensional Lie algebra with a basis {A m , Gi} me z i i e z >0 satis- 
fying the commutation relations: 

[A m , A n ] = 4G m - n , [A m , G{\ = 2A m _i - 2A m+ i , [G m , G{\ = , 

where G-i := — Gi and Go := 0. The elements Aq,A\ satisfy the Dolan-Grady (DG) relation [24]: 

[A u [A u [A^Ao]}} = 16[Ai, Aq] , [A , [A , [A , A ± }}} = 16[A , A x ] , (46) 

and the Onsager algebra is characterized as the Lie algebra generated by the DG pair {^o^i} 
[19, 38]. A useful realization of Onsager algebra to identify it with the Lie-subalgebra of the s/2-loop 
algebra, sfalz, z^ 1 }, fixed by a standard involution through the identification [38]: 

A m = 2z m e + + 2z~ m e- , G m = (z m - Z - m )h , m £ Z, 

where e^,h are s/2-generators with [e + ,e~] = h, [h, = ±2e ± . It is known that all finite- 
dimensional irreducible representations of s^lz, z~ l ] are given by tensoring a finite number of 
irreducible s^-representations through the evaluation of z at distinct non-zero complex values Oj's 
[16] . The Hermitian irreducible representations of Onsager algebra are obtained by passing through 
sfalz, z _1 ]-representations, with the further constraints on evaluated values a/s : 

\aj\ = 1, cij 7^ ±1, cij 7^ 1 for j / k. (47) 

By this, A m = 2J2] =1 (afe+ + aj m ej), G m = E"=i(a™ - a^ m )h r The a/s give rise to the 
reciprocal polynomial of degree 2n, 

2n n 

Y, a kZ k := H(z - aj )(z - aj 1 ), 

k=0 j=l 

characterized by the finite recurrence condition of A m 's: Efc=o a k A k+e = for t G Z [18, 23, 38]. 

It is known that the logarithmic derivative of transfer matrices T p (x q ,y g ) at the superintegrable 
p in (23) gives rise to the following Hamiltonian of Zjv-symmetric quantum chain in [32, 34] (see, 
e.g. [2, 3, 8, 37] and references therein), 

H(k') = H + k'H x (48) 
where Hq and H\ are Hermitian operators of <8> C defined by 

L N—l vn L N—l 711 7— n 

e=l n=l 1=1 n=l 

Note that for N = 2, this is the Ising quantum chain. The pair of operators 

A = -2JV -1 .Ho , A 1 = -2iV" 1 J fri 

satisfy the DG condition (46), hence gives rise to a Hermitian representation of Onsager algebra. 
By this, one can study the eigenvalue problem of H(k') in (48) through the representation theory 
of Onsager algebra. By [2, 8, 11] (or see the discussion in section 4 of this paper), only spin-i 
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^-representations occur in the associated Onsager algebra representations, hence H(k') has the 
Ising-like eigenvalues: 

n -. 

a + 0k' + N ±\j I + 2k' Cj + k' 2 , Cj = cos 9j = -(aj + aj 1 ) , a, G R . (49) 
For later use, we define 

AT— 1 vn AT— 1 vn y—n 

Ee= y^^, y e= yfl^±L : i < £ < l . (50) 

n=l n=l 

Then # = -2£ii Hi = -2E<=i 

Proposition 1 i^o commutes with T 2 (t) for all t G C. 

Proof. By .ffo = — 2£^ =1 H^, the commutativity of -ffo an d T 2 (t q ) follows from [S^,r 2 (t g )] = for 
1 < £ < L. By (25), it suffices to show the commutativity of Ep and (Gj)^(t) for all £,j,a,0. By 
the expressions of Eg in (50) and G^(i) in (24), one needs only to show Eg commutes with (G^)o(i) 
and (G/)?(t). By Z^JQ = and = 1, one has 

(l - wX t )Z t ~t = (J2n=o ^)(1 - = (s/ + E^o 1 (i - ^)^ 

= E e (l - uX e )Z e ; 

(i - xi)ZfE t = (El-o 1 - x e )zi l = (e £ - e^-o 1 A?) (1 " V 

= E e (l - Xf)Zj x . 
Therefore [(G,)J(t), E e ] = [(G,)?(i), S,] = 0. □ 



Proposition 2 i2i commutes with r 2 (t) for all t G C. 

Proof. By #i = -2^ the commutativity of #i and r 2 (t g ) follows from [V e ,T 2 (t q )] = for 
1 < £ < L. For a given £, commutes with (Gj)P(ut) for j / 1,1 + 1 in the definition (25) of r 2 (t). 
Here we use the identification L + 1 = 1. Denote by G^ + i(t) the 2-by-2 matrix with the entries, 

(G e ,e+it (a,/3 = 0,1), defined by (G w )£(t) := E j=0 ,i ((G<)&(G*+i)?) (t). Hence it suffices 

to show the commutativity of E^ and (G^ + i)^(t) for all a,/3. Using (24), one has the following 
expressions of (G^+i)f : 

(G e/+1 )° (t) = 1 - ^"^Iilf+i^^ + rj- 2 t 2 X e X e+1 

- V -H(x e (i - uZiZ^) + x e+1 (i - z.zjlS) , 

r/5(G v+ i)J(t) = Z i+1 -uj 2 XiXi +1 Zi-r]' 1 t(Zi-ujXiXi +1 Zi +1 ) 

- (rj-HXtil - uZeZg^) + uX e+1 (l - ZeZ^Sj Z e+1 , 
-vn-HGw+^it) = Zt 1 -u,X t X t+1 Z£ 1 -T ] -H(Z£ 1 -X e X e+1 Zf 1 ) 

- (xe(l - ujZgZ^) - Xe+i(l - ZeZ^^j Z^ 1 , 
(G e/+1 )\(t) = cj 2 X e X e+1 - r l -H(l+u J X e X e+1 )Zt 1 Z e+1 + V ~ 2 t 2 

- (xe(l — ujZgZ^) + u;Xe + i(l - ZtZg+ifj Z e 1 Zg + i. 



12 



By the expression of ^£ in (50), ^ commutes with Zi,Z^ + \ and XtXg + i. Furthermore, by 

* t X t {\ - uZeZ^) = X t (V t + E^=o w B ^)(l - ^^-\) = X t (l - uZtZ^jVt, 
V e X e+1 (l - ZtZ^J = X e+1 - En^o 1 ^+ n i)(l - = Xe+i (i - ^+1)^, 

all the terms appeared in the expressions of (G^+i)^ commute with Then follows r 2 (t q )] = 
0. □ 



Since Hq, H\ form the DG pair which generates the Onsager algebra representations, Proposition 
1 and 2 imply the Onsager algebra symmetry of T^ 2 \t). Note that the spin-shift operator X(= 
n£=i commutes with Hq and H±. By the expression (27) of we obtain the following result: 

Theorem 1 T/ie -operators in the superintegrable CPM possess the Onsager algebra symmetry, 
i.e., T^\t) 's commute with Hq and Hi, hence the corresponding Onsager irreducible representations 
give rise to the degenerated eigenvectors of r^\t) for 2 < j < N. Indeed, the eigenvectors of H(k') 
associated to the eigenvalues (49) form a degenerated eigenspace of r^\t) 's. 

□ 

Remark. By the commutative property (13) of and T p for the superintegrable element p, 
commutes with the Hamiltonian H(k') in (48). Since H(k') is the sum of Hq and k' H\ for a non-zero 
k' which is related to r/ in by rj = (j^p) 7r , Proposition 1 is equivalent to Proposition 2 when 
employing the theory of CPM. Here we give the algebraic verification of Onsager algebra symmetry 
of r^ 2 \t), hence r^-matrices, directly from the explicit form of t^ 2 ) -matrices and Onsager algebra 
generators, not making use of the commutative relation (13) of and T. However, to questions 
concerning the nature of the representations of Onsager algebra, one still needs to use the complete 
theory of superintegrable CPM, which we are going to discuss in the next section. □ 



4 Bethe Equation of r^-matrix, Degeneracy and Eigenvalues of 
r^- matrix 

In this section, we discuss in details the nature of Onsager algebra representations and r^)- 
eigenvalues in Theorem 1 using results essentially taken from works in [1, 2, 10, 11]. 

In the discussion of superintegrable CPM in section 1, it is natural to introduce the variable 

t = n~H 

where as before, rj = (y^pp) 1 ^- The hyperelliptic curve (6) becomes 

~ N (l-k'xm-k'x- 1 ) , , ,a + i, 2 t N -n- 2N 

f = ~ JT^) 2 or equivalents, (__)2 = __J__ , ( 5 1) 

i.e., ( ^pr5f~ ) 2 = ■ A PP!y the S au § e transform M- l G{t)M to G(t) in (24) with M = 

dia.[l, 772], then M~ 1 G(t)M has the expression 

~~ ( 1-tX {1-UX)Z\ 
° [t) ~ { -t(l-X)Z- 1 -t + coX ) ' 
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which is again a Yang-Baxter solution (8) with a = b = c = d=l. Then T^ 2 \t) in (25) can be 
written as T^ 2 \t) = T^(t), where 

L 

T {2 \t) = t W ((g) G e )(ut) , G e (t) = G{t) at Hh site. (52) 

£=1 

Write T^\t) = f^\t), then the fusion relation (26) becomes 

T^\t)T^ 2 \J- l t) = (1 - J- X t) 2L T { i-^(t) LO L X + 1 < j < N , 

5=^+1) (i) = (i_i)2L ^^-ij^^j + 2(1 - t N ) L . [bS) 

Note that the f ^-matrices are temperature ( i.e. k') independent. With operators Q cp , Q cp in 
(28), the relations (29)-(32) become the following r( 2 )Q cp , r^)Q cp -relations, 

f {2 \u~ l t q )q cp {q) = (1 - u- l t q ) L q cp {U- l q) + (1 - t q ) L Q> cp (Uq)u L X , (54) 

r {j) (t q ) = ]T ( n V° (1 ~^lf ^ Q C p(g)Qc P (^ m g)' 1 Q C p(^g)Q cP (^ m+1 g)' 1 (^^) J ' m ' 1 ), (55) 

m=0 V (1 - / 
and the Q cp Q cp -relation 

q^q^*) = ^ ( n^iW ^fo + nI:>-^)^ (JV " J)(a;JTg) ) ' - J - 7V; 

Qc P (g)Qcp(C9) = irfW(t 5 ) . (56) 

By [2, 10, 11], eigenvalues of the transfer matrix of CPM in the superintegrable case are solved 
by Bethe-ansatz method as follows. Let F(t) be a ^-polynomial (of degree m p ): 



F(t) (= F(t; Ufnj> )) = H(l + u Vi t) (57) 



where v i, . . . , t> mp are the parameters with {—Vi) N / 0, 1 and Ujt; - 1 / 1, w for i / j. Associated to 
F(i), we consider the following P-function 

,~ . p^- 1 (i - t N ) L (uj j ty Pa ~ Pb , , 

nt ) (=^ ;pi ,...,^))=^e ' (58) 

where P a ,Pb are integers satisfying 5 

0<P a + n<iV-l, P b -P a = Q + L (mod AT), (59) 

with Q the Z^r-charge as before. For arbitrary vi, . . . ,v mp , V(t) is a rational function invariant 
under the transformation: £ wt. The necessary and sufficient condition for V(t) to be a t- 
polynomial is the Bethe-ansatz-type constraint for parameters v\, . . . , v mp , ((4.4) in [2], (6.22) in 
[10]): 

i^^)^-*-*--^*^^, i = l,...,m p . (60) 

Vi + UJ 2 ^ Vi - UJVl 



'Here we consider only the case with periodic boundary condition by assuming r in [10, 11] equal to 0. 
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The above relation, up to the phase factor, is the Bethe equation of spin-^i XXZ chain for size 
L with the anisotropy 7 = and periodic boundary condition [10, 17]. An equivalent form to the 
above Bethe equation is 



(^^^ = -^-^11^^, i = l,...,m pj (61) 

Vi + UJ 1 fj^ Vi - UJVl 

and P a is the integer determined by (59). Then V{t) is a t^-polynomial of degree mg with zeros 
t N = si, . . . , s mE . Note that by the theory of the CPM V(l) / 0, indeed its value determines the 
total momentum P through the relation 



JP _ . .-p„ TT 1 + ^ _ ^(!) 



(see (2.24) of [2]). By e lPL = 1, the relation (62) yields the restriction of the non-negative integer 
m p in (61): 

LP fe = m p {Q- 2P b - m p ) (mod N), (63) 

(for iV = 3 case, see (C.3) (C.4) in [2]). Conversely, if m v satisfies the above relation, the right hand 
side of (62) is a Lth root of unity, hence determines the total momentum P. Furthermore the integer 
Pb in the theory is chosen to guarantee V(0) 7^ 0, (which will become clear later on). Therefore 
the t^-zeros of the polynomial V{t) have all Sj / 0, 1. For each j, let Xj be the complex number 

defined by the relation (51) with (t N ,X) = (sjjA^ 1 ), and denote ±uij = - = ( ^7^3 ) 2 - 
Hence Xj 7^ 1, /c /=tl , 0, 00, equivalently, Wj ^= 00, ±1. Define 

g(A) = n A + 1 *<*-'>»' ■ (64, 
i=i 

Then g(A; /=tl ),a(0) are non-zero. By \(X + 1 + (A - l)wj)(X~ 1 + 1 + (A" 1 - l)vjj) = one has 

the following relation between the P-polynomial (58) and Q{X): 

m = a(A)C?(A- 1 ). (65) 

By (1.11) in [2] and (21) in [11], one has the following T(g)-eigenvalue 6 : 

T (i) = n l £Z%fX )L (v^ q ) p ^y q f^ p ^G(X q ), (66) 

and T(q) is equal to e lP T(q) with the expression 



6 Here we modify the expression (21) of [11] by a multiple factor e~ lP / 2 so that T(q) takes the value 1 for q — p, 
which would be more appropriate for the consistency with the analysis made in formula (1.11) in [2], elucidated later 
on in this section. Note that the expression here agrees with the conjectured formula (2.22) of [2], in which the 
notations Af, A, 77, NP C , wt correspond respectively to L, k' , r/~ , P M , 1 ~ fc we used in this article. Indeed the T q in 

formula (21) of [11] is equal to N L £J,^%l (^^/'(^y/V^ ^ ^l^ iK)- 
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Hence we have the following eigenvalues for operators Q cp (q) and Q cp {q)' 
Qc P (q) =(v^x q ) p «(r ) ^y q ) p ^- P ^y(\ q ), 

which imply Q cp (p) = 1, Qc P (p) = e lP l. As q tends to p, by setting x q = 772(1 — 2fc'e + 0(e 2 )) with 
small £, to the first order one has y q = 772 (1 + 2fc'e) , /x 9 = 1 + 2(A/ — l)e, and formula (1.11) in [2] 
for the CPM-transformation matrix 

T(x q , y q ) = 1[1 + e(N — 1)L] + eH{k') + 0(e 2 ) 

where is the Hamiltonian (48). Then the T(g)-eigenvalue (66) gives rise to the following 

energy value of H(k') ((2.23) of [2]): 

/ \ rn E 

E = 2P^ + Nm E - (TV - 1)L + k'({N- 1)L - 2P^ - Nm E + 2(P b - P a ) ) + iV(l - k') ]T ±u>;- (69) 

Write -ff(fc') = ^-(Aq + fcM.i) as before with Hermitian ((g) C Ar )-operators Ao, A\ satisfying the 
DG relation (46), hence generating a Hermitian representation of Onsager algebra. Each Onsager 

L N 

algebra irreducible subrepresentation of g> C is obtained by an evaluation representation of 
sl2[z,z~ 1 } by evaluating z at finite complex numbers a,fs satisfying (47), then composed with 
irreducible representations of the ^-factors, characterized by the spin aj € ^Z>o- Through this, 
one can write A = =^-a + 2X^(e+ + ej), A 1 = =Jff3 + 2^(oje+ + a^ej), which imply the 
eigenvalues of H(k') with the form 



a + (3k' + 2NJ2 mj^l + 2k> cos^ + fe' 2 



where a,/3 G R, cos6*j = a,j + aj 1 , and mj runs through —aj,—aj + l,...,<jj [18, 19, 23, 38]. 
Comparing the above formula with the energy expression (69), we find 

1 1 + s ■ 

aj = -, cos9j = - for all j , (70) 

2 1 — Sj 

then follows the expression (49). Since cost's are distinct real numbers with |cos#j| < 1, the 
r^-polynomial associated to V(t) in (58) is simple with negative real roots Sj (1 < j < mg). In 
particular, 'P(O) 7^ 0. Hence we have shown (i) of the following theorem except the formula of mg. 

Theorem 2 (i) Let {i^}™^ be a solution of the equation (61) for L,Q,Pb with the relation (63), 
and V(t) the polynomial (58) associated to {vi}™? 1 . Then V(0) 7^ 0, and Vit) is a simple t N - 
polynomial of degree tue with negative real roots. Furthermore, wie = [ — ^ L ^ Pb 2mp ] where 
[r] denotes the integral part of a real number r. 

(ii) For {vi} r ^F 1 in (i) with the total momentum P, denote v[ = a; -3 ^" 1 , and Q' = —L — Q, 
Pi = —Pf, — m p (mod N). Then {v! i }™T 1 is a Bethe solution of (61) for L,Q',Pl with the total 
momentum P' = —P. Furthermore, the polynomial V(t;v[, . . . ,v' mp ) in (58) associated to {v' i } 1 ^ 1 
is related to V(t; v±, . . . , v„ lp ) by the following reciprocal relation, 

V(t; v[,..., v' mp )u,n l[ v > = F m * Pit- 1 - Vl , . . . , v mp )^ J] ^ • 
i=i i=i 
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Proof. It is easy to see that the < [ ( JV ~ 1 ) L ~-^~~ Pi >~ 2m p ] - L e t Q' ; p^ p> be those defined in (ii) in 
connection with i> j, Q, P&, P in (i). One can easily see that {-Uj}™^ as a Bethe solution of (61) lq^ 
is equivalent to that for {t^}™^ in (61)x, q' p'j the same for (P, {^j}^) and (-P',{^}^fi) with the 
relation (62). By (59), P b -P a = Q+L, P b '-p] = Q'+L (mod iV), and < P a +P b , P^+P b f < N-l. 
Hence P' a + P b ' = —L — P a — P b — 2m p (mod N), which implies 

r L-P a -P b ~2m p ] _ r L-P' a -P^-2m p ] _ -L-P a -P h -2m v -P' a -P' h 
L N J L N J AT ' 

and [ ( iV-l)L^-P b - 2mp ] = (AT-llL-P.-^-^-i*-^ ug . ng rdation 



(1 _ v > _ ^-l)L- 2mp -P a -P 6 -P^P^ 1 _ f-AT^+P, jj** 



(i - o L ((i + ^)(i + u 2 v'S) (i - nr= p i ((i + ^-^(i + ^t- 1 : 

one has 

V(t;v[,.. .,v' mp )uj p b \\v[ = PI iv 1 V(t- l -v u .. . lVmp )uj Pb J] Vi . 

i=l i=l 

Set t = in the above relation, then V(0; v[ , . . . , v' m )w p ' n™ P i ^ = Cuj Pb ]~I^=l v u where C 

is the coefficient of r 1 jv J in the polynomial V(t; v±, . . . , v mp ). By (i) for V(0) 

(associated to v^s) , we have V(0; v[, . . . , v' m ) / 0, hence C / 0. Then follow the equality 
rriE = [ ( N ~ 1 } L ~ p ^~ p b- 2m p ] i n anci the reciprocal relation in (ii). □ 

Remark. The dominant eigenvalues of T(q) are in the case P& = m p = [11], where tue (= 
[ ^ Ar ~ 1 ^ L ~ Pa ]) = [ ^ JV ~^- L ~ < ^ ]. The mathematical structures of polynomials P(t) and the correspond- 
ing polynomials in Onsager algebra representations were discussed in [30, 33, 40] □ 

Let {vi}™T 1 be a Bethe solution of (61), and F(t) the polynomial (57) defined by Vj's. Associated 
to F(t), one has the Q cp and Q cp -eigenvalues (68), with the corresponding r^-eigenvalue for all j. 
Then relations (54) and (59) yield ((6.18) in [10]): 

TV>{u>-%)F(t q ) =u~ p *(l -u- l t q ) L F(u~ l t q ) + u p »(l -t q ) L F(ui q ) . 

Set t q = —(uJVi)" 1 in the above relation. One obtains the Bethe equation (60), equivalent to (61), 
of ViS with the following polynomial expression for f (^-eigenvalue: 

T (t)(-T (t, Vl ,...,V mp ))- ^UZ\(l+^ Vl t) 

(71) 

By (55), one has the r^-polynomial expression: 

Then the following relations ((6.17) and (6.24) in [10]) 7 hold: 

F°+ p »F(t)F(uH)V(t) = "7? {1 - TN) ~ r fU)(t) + "T: {1 - JN) ~M N -i\un), 0<j<N, 

?W(t) = t p «+ Pb F(t) 2 V(t) (for j = AT) 



7 The formulae here differ with those in [10] by some scalar factors as we use conventions in [2, 11] in this article. 
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where V(t) is in (62). Note that by (65), the above relation is equivalent to the Q cp Q cp -relation 
(56). Hence we have shown the following results: 

Theorem 3 (i) The Bethe equation of t) for the L-size t 1 ■ '-matrix (52) in the six-vertex 

model (with external fields), is the relation (61) (or equivalent to (60),), where the equation depends 
on two quantum numbers Q,Pb £ T*n with the constraint (63) on non-negative integer m p . The 
total momentum P is given by (62). 

(ii) Let {vi}^ be a Bethe solution of (61), and F(t)(= F(t; v±, . . . , v mp )) the polynomial de- 
fined in (57). Then there associates an eig en-polynomial T^\t) (= f^(t; v\, . . . , v mp )) of the f^- 
operator for 2 < j < N, given by (72) such that the relation (73) holds. 

□ 

Remark. One may as well obtain polynomials r^)(t)'s by using the (t)-expression (71) and 
the fusion relation (53). □ 

Note that Theorem 3 was derived from Q cp -eigenvalues in (68), where all possible Q(X q ys give 
the same f ^-polynomial for 2 < j < N. By (64) and (69), the eigen- vectors corresponding to those 
available Q cp -eigenvalues form a 2 mE -dimensional vector space which is the representation space for 
an irreducible Onsager algebra representation with the evaluation parameters determined by the V- 
polynomial (58) via the connection (70). Hence the Onsager algebra describes the full symmetry of 
the degeneracy space for an given ^-polynomial. This can be regarded as an explicit formulation 
of the Onsager algebra symmetry of r^-matrix in Theorem 1. As a corollary of Theorem 2 (i), we 
have the following result. 

Proposition 3 Let {fj}™^ be a Bethe solution of (61), f^\t) (= f^\t); v\, . . . , v mp ) the eigen- 
polynomial in (72) for 2 < j < N. There is an irreducible Onsager algebra subrepresentation space 

(of® C ) of dimension 2 l n > with T^'yt) as the t (j > -eigenvalue. 

□ 



5 Concluding Remarks 

In this paper we have shown the Onsager algebra symmetry of t( j ) -matrices in superintegrable 
CPM by using the explicit form of r^ 2 ) -matrices and Onsager algebra generators. Furthermore, by 
employing the theory of CPM, the nature of Onsager algebra representations and the full symmetry 
structure are revealed in the Bethe Ansatz approach using the explicit eigenvalues of superinte- 
grable CPM. Note that in this paper, for simplicity, we have confined our study of symmetry 
in superintegrable CPM only in the case of periodic boundary condition. We know that general 
skewed boundary conditions are required for the natural completion of the eigenvalue spectrum 
even in the Ising case. However, by results in [10], the discussions of superintegrable CPM in this 
article can be carried over as well to skewed boundary conditions in general, hence the Onsager 
algebra symmetry would be valid for the complete spectrum of r^-model. One may view the 
Onsager algebra symmetry of r^-matrix in CPM as the counter part of the s/2-l°°P algebra sym- 
metry of the root-of-unity six- vertex model found in [20] within the scheme proposed in [29] on the 
(conjectured) analogous theories between the eight vertex model and CPM. In this work, we have 
performed a rigorous study of the "P-polynomial (58), the pivotal ingredient in r^-degeneracy, by 
explicit results of eigenvalues of CPM transfer matrix. In the context of s^-loop algebra symmetry 
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of six- vertex model [22, 27], one has the Drinfeld polynomial, which plays the similar role of V- 
polynomial in CPM, however is not fully understood till now. The exact results we have obtained in 
this work could be valuable to pave the way to solution of some conjectured symmetry problems in 
the root-of-unity six-vertex model. A process along this line is now under consideration. Moreover, 
we have made a detailed analysis on the comparison on fusion matrices and functional relations 
of the root-of-unity eight-vertex model and superintegrable CPM. As a consequence, we found the 
relation (40) in eight-vertex model corresponding to TT-relations (31) in CPM, which has played 
the vital role for the derivation of all functional relations in [14]. However the quest for rigorous 
justification and physical interpretation of (40) has been left a challenge for the root-of-unity eight- 
vertex model as a completely parallel theory to CPM. The use for deeper understanding of the 
Q-operators would be expected. Nevertheless, the striking similarity of these two theories would 
suggest the universal role of CPM in regard to the "(^"-operator in Baxter's TQ-relation, which 
appears also in other solvable lattice models. We made a attempt in Sec. 2 to explore the common 
features of functional equations for other plausible, but hypothetical, lattice models, based on the 
structure of superintegrable CPM listed in Sec. 1. The principle could be very sketchy, and in 
part just based on speculation; however, we hope it might help to give a feeling of how the current 
research on symmetry of r^-models is oriented. 
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